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Abstract. - We investigate coherent backscattering of light by two harmonically trapped 
àtoms in the light of quantitative quantum duality. Including recoil and Doppler shift close to 
an optical resonance, we calculate the interference visibility as well as the amount of which-path 
information, both for zero and finite temperature. 



Introduction. - A wealth of information about the motion of microscopic partides can 
be gathered by scattering a well-controlled probe, typically light in the case of àtoms with 
optical resonances. Young's històric double slit experiment demonstrated that the probe light 
wave shows interference if it can propagate along more than a single path. Interference fringes 
therefore can give sensitive information about the scatterers, as was shown by Eichmann and 
coworkers [1] by realizing Young's double slits with two trapped àtoms. On the other hand, 
the motional state of the probed particle gives information about the probe field, as recently 
discussed by Eschner [2]. As a commonly accepted rule, one can have either full interference 
contrast or full which-way-information at the same time. Between these two extremes, there 
are interesting intermediate situations that can be quantified using general quantum duality 
relations derived by Englert [3] that have been used for instance in atom interferometers [4]. 
In mesoscopic samples of wcakly disordered clouds of cold àtoms, coherent múltiple scattering 
of light leads to enhanced backscattering [5,6]. 

In this Letter, we investigate coherent backscattering (CBS) of light by two àtoms that 
are trapped in harmònic oscillators. This model system can describe how atòmic motion 
destroys múltiple scattering interference via recoil and Doppler effects. We calculate the CBS 
interference visibility both at zero and finite temperature in shallow traps that allow to treat 
the limiting case of freely moving àtoms. The CBS double scattering geometry realizes a 
two-way interferometer where the two mutually exclusive alternatives are the order in which 
the photon visits both àtoms. Which-way information is then present if one knows which of 
both àtoms has been visited first. We calculate the which-way distinguishability proposed by 
Englert and show how it is physically encoded. 
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Fig. 1 - (a) Two harmonically trapped àtoms backscatter a probe photon along way A or B. (b) Recoil 
mechanism: which-way information can be extracted if the wave-packet displacement Ax — v lac At of 
the first scatterer during the total scattering time At is larger than the initial position uncertainty A. 



The setting. Consider two àtoms trapped in identical harmònic oscillators at fixed 
positions well separated from each other by many wavelengths of the probe light. A single 
incident photon with wave vector k m and polarization èi n is then scattered by the àtoms and 
detcctcd in the backscattering direction fe out = —k- m with polarization analysis in the helicity- 
preserving channel è ou t = è* n . A photon scattered in the backscattering direction by a single 
atom with a non-degenerate ground state has the same polarization, but opposite helicity, 
and does not contribute to the detected intensity. Since the àtoms are far from each other, 
the probability for repeated scattering is very small. The detector then receives a photon that 
has been scattered exactly once by each atom, either along way A or along way B, see fig. ^ 
The total amplitude is the coherent superposition of both amplitudes. Such an interference of 
counter-propagating múltiple scattering amplitudes in the backscattering direction is known 
as coherent backscattering (CBS) [5]. 

Without interaction, the free electromagnètic field and the infernal and external (motional) 
atòmic degrees of freedom are described by the Hamiltonian Hq = ií p h + -ffint + -ffext or 

H = H^I^ + huJ (Pel + Pc 2 ) + Muo^l + N 2 ). (f ) 

The first term constitutes the Standard Hamiltonian for free photons [7] . P C i = I ^ rn i)(^ m i I 
is the projector onto the degenerate multiplet J c = 1 of excited states with transition fre- 
queney ujq from the non-degenerate ground state J g = of àtoms i = 1,2. The atòmic motion 
is described by the number operators Ni = Ni X + Ni y + Ni Z = aja, of excitations in the 
three-dimensional isotropic harmònic oscillators with frequeney Lüho- 

The dipole interaction V = —D\ ■ E(R±) — D 2 ■ E{Ri) couples the photon field both to the 
electrònic states and the motional degrees of freedom. Di is the electrònic dipole transition 
operator for atom i. The elèctric field operator E(R) is evaluated at the atòmic center-of- 
mass position Ri = Rf 1 + U{. The displacement operator Ui — \\ lo {a\ + o») from the trap's 
origin R^ measures distance in units of the oscillator length Aho = V h/ (2mL0h o ). The vector 
joining the two àtoms will be denoted R21 — R2 ~ Ri, with a distance Rq = R^ — R^ 
between the traps such that fem-Ro 3* 1- 

Transition operator for photon double scattering. - Since the àtoms are well separated, 
the total double scattering process is described by the product of two individual scattering 
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processes with free propagation inbetween. In the far-field photon propagator e lk I - ** 21 ' /I-R21 1 , 
we expand the absolute distance to linear order in the small relative displacement: I-R21I ~ 
i?o + Ro ( u 2 — «i). In the exponential, the zeroth-order contribution Rq drops out from 
all interference quantities, whereas the linear term generates a phase difference between the 
amplitudes of A and B and must be kept. To this leading order, the denominator can be 
takcn constant. 

The transition operator [8] for way A takes a useful form in time representation and 
interaction picture where Rj(t) = e tHat / h Rje~ lHat / h . Up to irrelevant prefactors, 

poc pOC 

T A = / e il(s+t) e -ik oM R2(0) e ik'R 2 (-s) e -ik'R 1 (-s) e ik i nRi(-s-t) ^dg. (2) 

Jo Jo 

Read from right to left, it describes how the photon is scattered first by atom 1, then prop- 
agates with fe' = fcm-Ro, and is fmally scattered by atom 2. The transition operator Tb of 
the reverse way B is obtained from Ta by the substitution Ri <-> i?2,fe' — > —fe'. The com- 
plex detuning 7 = Lo- m — ujq + iT/2 of the probe frequency Wi n from the transition frequency 
includes the spontaneous decay rate or inverse lifetime T. Retardation times of order Rq/c 
have been neglected in the time arguments of the operators. Indeed, the photon scattcring 
by two àtoms defines two distinct time scales for the atòmic motion: first, the total inverse 
detuning | ^ | 1 which is simply 2/T at resonance. Second, the free propagation time Rq/c 
from one atom to the othcr. For resonant atòmic scatterers with 1/T of order 10 7 s or larger, 
and typical distances Rq of lmm or less, one has Rq/c <C 1/T such that the scattering is 
resonance-dominated, and the frec propagation time can be safcly ncglcctcd. In all of the 
following, we have in mind the limit of quasi-free àtoms and therefore consider the case of 
shallow traps uj\ io <C T in which the oscillation period of an atom is much larger than the time 
r _1 it takes to scatter a photon. 

Quantitative quantum duality. - The incident photon can choose between two a priori 
indistinguishablc paths, way A or B. Following Englert's fashionable choice [3], we call this 
binary degree of freedom a qubit. Let | A) and |B) denote the choice of way A and B. Scattering 
entangles the qubit with the motional degrees of freedom. Consequently, the atòmic oscillator 
states can serve as a which-path detector. 

Prior to scattering, the total initial state of qubit and detector is p^t = p AB ® pW. 
The qubit is in the pure state p^ B = | (|A) + |B)) ((A| + (B|), the symmetric superposition of 
cqually probable ways. An external cooling làser field serves as a thermal bath for the trapped 
àtoms. The detector is thus in a thermal state pW — e~ l3HaM /Z at inverse temperature 
fi = 1/fceT with partition function Z = tr{e~^ ffcxt }. The total final state is then obtained by 
applying the transition operators associated with way A and B: 

(& = T^T, T=\A)(A\®T A + \B)(B\®T B . (3) 

In general, T is not a unitary operator since it describes only the scattering amplitude around 
the backscattering direction. The factor Ia + Ib = tv{TAp^T A } + tr{TBP^T^} guarantees, 
however, that pl ot is properly normalizcd. 

Adapting Englert's general definitions [3], we can express the visibility V and the distin- 
guishability V obeying the fundamental duality relation V 2 + V 2 < 1 as follows: 



V = 2 



V = 



tr{T B/ 9«T{} tr TaP^TÍ - T B p®T} 



Ia + Ib ' Ia + Ib 



(4) 
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The visibility is the ratio of interference contribution and background intensity Ja + and 
therefore equal to the CBS contrast. The distinguishability T> describes the maximum which- 
way information available in principle, i.e., that can be extracted by the additional measure- 
ment of an optimal detector observable; tr|X| = tr{V X^X} denotes the trace-class norm of 
the operator X. Another interesting quantity is the predictability 



V 



tr{T A p«Tl-T B p«Tt}| |Wfi| 

(5) 



Ia +Ib Ia + Ib 



The predictability measures the amount of which-way information available a priori, as for 
instance in unbalanced interferometers like Young's double slits with different widths. For 
balanced interferometers Ia = Ib, the predictability vanishes, V = 0. 

With the hclp of these quantities, we can quantify the breakdown of coherent photon 
backscattering by mobile àtoms, both at zero and finite temperature. The distinguishability 
T> is difficult to evaluate in our case because it is defined via the trace-class norm of operators 
on the infinite-dimensional Hilbert space of harmònic oscillators. On the contrary, visibility 
V and predictability V involve a thermal harmònic average of products like t[t b or T A T A . 
Since the transition operator contains only exponentials linear in the displacement, taking 
the trace amounts to Gaussian integration. 

Zero temperature. At zero temperature, the àtoms are initially in their respective 
harmònic oscillator ground states. In other words, the detectors are prepared in pure states. 
Setting = in J3J permits to show that in this case the duality relation is always 

saturated: V 2 = 1 - V 2 . Thus, we can use the visibility V, much easier to calculate than the 
distinguishability 2?, in order to understand how much which-way information is present and 
how it is encoded. Carrying out the thermal average, we find in the limit T — > 

v 2 = eVm) + 2C 2 K0 = i-v 2 , (6) 

V 2 = g 2 (uj in h (7) 

neglecting higher-order terms 0(Ç i ,x 3 ,( 2 x) m the two relevant small parameters defined as 
follows: The influence of atòmic recoil without harmònic trapping is encoded in the factor 

with ?ïa; rec = h 2 k 2 1 /(2m) the recoil energy, S = uj in — lu the probe detuning from the atòmic 
resonance, and 7 = S + iT/2 as before. The harmònic trap enters through the parameter 

C M = | 7[ 2 ■ ( 9 ) 

It is worthwile to discuss the physical significancc of both parameters. 

Let us first consider a very shallow trap cjho *C w rec where the free recoil effect described 
by g(u>i n ) dominates. At least to order \ 2 , w e find V 2 + V 2 = 1. Together with the general 
property V < T> (read: the a priori information cannot be larger than the total available 
information) this implies that all which-way information is actually available a priori: T> = V . 
Just as in the case of other asymmetric interferometers, this predictability is due to unbalanced 
scattering amplitudes. At perpendicular scattering i? ■ fein = 0, the situation is completely 
symmetric, and both paths are equally probable. But for the extreme case Rq ■ k[ n = ±1 of 
àtoms in line with the probe, the atom in front scatters either the incident probe photon at the 
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Fig. 2 - (a) Frequency-shift mechanism: a small frequency change Aui on the steep slope of a resonance 
curve implies a large scattering cross section difference Aa that can provide which-path information. 
(b) CBS visibility @ as function of the thermal Doppler shift I1U1 at resonance 5 = 0. 



làser frequency or the already scattered photon on its way out again, now at làser frequency 
minus twice the recoil. A different frequency generally implies a different scattering cross 
section such that the two paths have different probabilities. For a given position configuration, 
this information is known a priori without the necessity to perform any measurement. More 
quantitatively, the relative change in the resonant cross section a(S) = <ro[l + (25 /r) 2 ] -1 under 
a small frequency change Auj -C T thcn is | Acr/ cr| = Auj2S/\j\ 2 , see fig. |2ta). For the actual 
frequency change Aw = 2(Rq ■ k- m )u! rcc , one finds exactly |Aer/er| = g>(c<Ji n ). To this order, the 
predictability vanishes at exact detuning, since a small frequency change on the flat top of 
the resonance Lorentzian has no effect. 

Lct us now intèrpret the influence of harmònic trapping at zero temperature. The pre- 
dictability J7J) contains no contribution in £ 2 (íJi n ), which indicates that this parameter encodes 
which-path information that may only be revealed a posteriori by an appropriate measure- 
ment on the detector. In a temporal picture (see fig. ^b)), one can determine which way 
the photon has taken if one can measure which atom has scattered the photon first. This is 
only possible if the initial position uncertainty A is smaller than the distance Ax travelled by 
the first scattering atom. More quantitatively, the first scatterer takes up the recoil in the 
direction of k ul and travels the distance Ax = v ICC At — \/2huj ICC /mAt during the whole time 
Aí = 2 1 | 1 until the emission of the final photon by the second atom. ( 1 ) At that moment, 
the second atom receives the same recoil in the direction fei n , such that the final momentum 
states for both ways become indistinguishable (note that the recoil in the direction joining 
the àtoms is exchanged instantaneously, because we can neglect the propagation time). The 
zero-temperature position uncertainty is the oscillator lcngth Aho, which finally shows that 
which-path information is indeed present on the scale Ax/X^o = 2Ç(u>i n ). 

Finite temperatures. - If the àtoms are coupled to a thermal bath, the detector is not 
prepared in a pure state, and one expects that the duality relation no longer saturates. We 
thus have to calculatc visibility and distinguishability separately. 

Aftcr the Gaussian thermal average, the visibility is given in terms of a product of am- 
plitudes like for way A and B. It can be evaluated numerically for all vàlues of the de- 
tuning S and of the thermal trap parameter £ 2 (T, Ui n ) = coth(i/?Stt'ho)C 2 ( w in)· This param- 

( 1 )At is simply twice the total contribution per atom that can be justified by a stationary-phase argument 
if both the Wigner time delay of the phase and the amplitude change of the resonant scattering function 
t(<5) = íq/(<5 + *r/2) are taken into account. 
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eter corresponds to the zero-temperature trap parameter (JÜJ) at effective oscillator excitation 
[2n(T) + l]who, with n(T) = [ e ' 3riüJh ° - the Bose-Einstein distribution function. For highly 
excited àtoms Phw^o <C 1 the oscillator frequency drops out of the defmition of the thermal 
trap parameter which becomes 

U = ^j^, (io) 

the classical thermal Doppler shift (v 2 ms = (v 2 )) in units of the linewidth. Figure^b) shows 
a plot of the visibility at exact resonance S = as function of ^ on a semi-logarithmic scale. 
Hcre, we concentrate on the effect of harmònic trapping at finite temperature and neglect the 
free recoil contribution, which is justificd if x/£ 2 j ^ 1- The visibility decreases monotonically 
with temperature. For ( c 2 ] < 1, we can expand to lowest order and find analytically (compare 
with eq. ©): 

V 2 - 1 - 2£ 2 , (11) 

At this point, we recover the case of free thermal àtoms. Our result (jl·ljl agrees with the CBS 
contrast calculated in [9] for the low-temperature case £ c i 1 (if the avarage medium effect 
included there is disregarded). Note that our calculation is also vàlid in the high-temperature 
regime £ c i ^> 1 where the visibility goes to zero as V ~ £~ 2 . 

It is too difficult to evaluate analytically the trace-class norm for the distinguishability 
D with the transition operators in the general form (J2J). Therefore, we formally expand the 
exponentials in powers of the Lamb-Dicke parameter &i n Aho = V uj icc /üj\ w . The leading order 
expression for T> in the shallow-trap limit [3huJho *C 1 becomes 



« 2V2|í| 

\l\ 



(12) 



Here, p — (h/2i\\ lo )[à — à^] is the momentum associated with à = (ai — a2)fci n , the anti- 
symmetric oscillator mode projected onto the probe direction fei n . It is reasonable that only 
this mode should be relevant: symmetric motion cannot encode differential information about 
the paths (and could be disposed of by transformation into a co-moving frame), whereas the 
momentum along the only other available direction Rq is exchanged instantaneously. To lin- 
ear order in = fíhui-^o, the positive operator V = e _/3a a can be moved outside the absolute 
value (using \Vp\ ~ y (V^V)(p^p) w \V\ \p\ since p and V commute to zeroth order in (3). This 
makes the distinguishability proportional to the thermal expectation value of the momentum 
modulus, which is easy to evaluate for a Maxwell-Boltzmann distribution: 

/ -\S\2v rec (\p\) 2 \S\ 
^ = v2— — — = —= — Ç cl (U) 

It can be read as the relative shift T> = 7r _1 / 2 |Acr / 'a\ of the resonant cross section under the 
small Doppler frequency change Ato — 2ki n v Tms . This suggests an interpretation in terms 
of the frequency-shift mechanism (cf. figEt a )) : measuring the actual atòmic velocities 
a posteriori, one can determine which way was more probable: if atom 1 is found to move 
towards a red-detuned probe and atom 2 to move away from it, way A was more probable 
than B because atom 1 had a larger absorption cross section — this is simply the principle of 
Doppler cooling. Contrary to the zero-temperature case, now the predictability V is zero to 
the order considered because the interferometer symmetry is reestablished by the isotropic 
thermal average. Naturally, as in the case of harmònic trapping at zero temperature, the 
distinguishability vanishes also at exact resonance (5 = because no differential information 
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can be encoded on the flat top of the resonance Lorentzian. In the high-temperature limit 
£ c i » 1, i.e., for very large frequence shifts Alü 3> T with moderate detuning, one explores 
the flat wings of the Lorentzian, and the distinguishability must then drop to zero. 

At finite temperature, the quantum duality is no longer saturated and reads to order 



This value of V assumes an optimal measurement within just the system of center-of-mass 
motion. But the initial thermal distribution rcquires the presence of a thermal bath (e.g., 
cooling làsers) . The total system (detector plus bath) is described by a pure state for which the 
quantum duality saturates. However, the which-path information could only be retrieved by 
an optimal measurement including all bath degrees of freedom that are not undcr our control. 
Measurements on just the detector necessarily correspond to a non-optimal measurement with 
respect to the whole system and thercfore imply a reduced distinguishability. 

Summary. - We have studied the coherent backscattering effect from two trapped àtoms, 
an interesting case study for quantitative quantum duality in a physically realistic setting. At 
zero temperature, inelastic scattering due to the recoil effect provides a priori which-path 
information. Furthcr which-path information can be gained by measuring which atom has 
scattered the photon first, which is possible if the initial position uncertainty is small enough. 
At finite temperature, the atòmic movement destroys the interference once the average Dopplcr 
shift becomes of the order of the resonance width. Which-path information resides in the 
atòmic velocities via resonance conditions, but the duality inequality no longer saturates 
because an initial thermal state corresponds to a non-optimal detector preparation. 

We hope that these considerations may stimulate further work, both experimental and 
theoretical, on the interference probing of quantum dynamics of trapped partides. 
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